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^ '■ 

, We investigate the rectification of an ac bias in Luttinger liquids in the presence of an asym- 

(•^^ ■ metric potential (the ratchet effect). We show that strong repulsive electron interaction enhances 

' the ratchet current in comparison with Fermi liquid systems, and the dc I — V curve is strongly 

I > , asymmetric in the low- voltage regime even for a weak asymmetric potential. At higher voltages the 

■ ratchet current exhibits an oscillatory voltage dependence. 

O: 

, Asymmetric conductors have asymmetric I — V curves. This phenomenon is known as the diode or ratchet effect 

' and plays a major role in electronics. Recently much interest has been attracted by transport asymmetries in single- 

i— ^ ' molecule devices and other mesoscopic systems 0. The idea that asymmetric molecules can be used as rectifiers is 
rather old 2], however, it was implemented experimentally Q only recently. Another experimental realization of a 
i-S^ , mesoscopic rectifier is an asymmetric electron waveguide constructed within the inversion layer of a semiconductor 
Jrj ' heterostructure . The ratchet effect was observed in carbon nanotubes ^nd strongly asymmetric I — V curves 
O , were recently reported for the tunneling in the quantum Hall edge states These experimental advances have 
^ ' stimulated much theoretical activity 0, S ^ liol Hlj with the main focus on the simplest Fermi-liquid systems |0| . 
^ , Transport in one-channel quantum wires, where electrons form a Luttinger liquid, differs significantly from the 
^ Fermi liquid case. In particular, impurity effects are stronger in Luttinger liquids, and even a weak impurity potential 
^ may render the linear conductance zero at low temperatures 13] . In this Letter we investigate the ratchet effect in 
^ Luttinger liquids. We show that strong repulsive electron interaction enhances the ratchet current, and the low-voltage 
(— ( part of the I — V curve is strongly asymmetric even in quantum wires with weak asymmetric potentials. 
Q ' We consider the ratchet effect in the presence of a weak asymmetric potential U{x) ^ Ep, where Ep is the band 
O i width. We calculate the ratchet current Ir{V) = [I{V) + I{—V)]/2 for a one-channel quantum wire with spin-polarized 
electrons. The ratchet current vanishes for systems with symmetric I—V curves. It can be measured as the dc response 
to a low-frequency square voltage wave of amplitude V. First, we consider voltages V < Vq = fivp / (ea), where vp is the 
^ Fermi velocity, e the electron charge, and a the size of the region containing the asymmetric potential. We find a weak 
0\ ratchet effect in the interval eVo > eV > y/UEp for both Fermi and Luttinger hquids, Ir ^ {e/h)U'^{eV)'^3 / E'p^+'^ , 
where g — 1 for Fermi liquids and g < 1 for Luttinger liquids with repulsive interaction. However, at strong repulsive 
^ interaction (the Luttinger liquid parameter g <C 1) and sufficiently low voltages, the ratchet current Ir{V) grows as 
^— I the voltage decreases until Ir{V) becomes comparable with the total current I{V) at eV = eV* ^ {UEp^y^^^^s\ 
At Ep ^ eV > eVo the ratchet current oscillates as a function of the voltage and can become comparable with the 
total current I{V) for any repulsive interaction strength. We also briefly discuss the ratchet effect in the presence of a 
strong asymmetric potential U > Ep. The complicated ratchet-current behavior is caused by the energy dependence 
of the effective impurity strength in Luttinger liquids This introduces an additional energy scale V* absent in 
Fermi-liquid systems. 

■ One-channel quantum wires can be described by the Tomonaga-Luttinger model with the Hamiltonian 

a : 

O . I dxl - nvp[4,l{x)idMx) - i'lix)id^Mx)] + U{x)p{x) + I dyK{x - y)p{x)p{y) \, (1) 



o 

where f/'/j and i^j^ are the creation operators for right- and left-moving electrons, ■0^ = V"]? + V'l gives the conventional 
' electron creation operator, p = ■0^';/' is the electron density, U{x) is the asymmetric potential, and K(x — y) the 
^ , interaction strength. Our aim is to calculate the current J as a function of the applied voltage V. We assume that 
- - ' the long-range Coulomb interaction is screened by the gates so that K{x — y) decreases rapidly for large [x — y). 
Electric fields of external charges are also assumed to be screened. Thus, the applied voltage reveals itself only as the 
difference of the electrochemical potentials El and Epi of the particles injected from the left and right reservoirs. 

We assume that one lead is connected to the ground so that its electrochemical potential Eb. = Ep is fixed. The 
electrochemical potential of the second lead El = Ep + eV is controlled by the voltage source. In such situation a sym- 
metric potential U {x) is sufficient for rectification. For example, in a non- interacting system I{V) ~ J^^ [1 - R{E)]dE, 
where R{E) is the reflection coefficient. If the only relevant scale for the energy dependence of the refiection probabil- 
ity is the band width - Ep then the ratchet current Ir ~ J^^ dE[R{Ep - E) - R{Ep + E)] w -2 ^ dER'{Ep)E - 
R{Ep){eV)'^ / Ep ^ U^{eVy/Ep for small U and V, and any coordinate dependence U{x). 
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A 'not-trivial' ratchet effect can be observed when the injected charge density is voltage- independent, E]^/^ = 
Ep ± eV/2. Symmetry considerations require an asymmetric U{x) for a non-vanishing ratchet current in this case. 
Also an electron interaction must be present. Indeed, for free particles the reflection coefficient R{E) is independent 
of the electron propagation direction jlj| and hence I{V) = —I{—V). 

The 'non-trivial' ratchet effect is absent in the first two orders in U{x). Indeed, in the lowest two orders the 
ratchet current 1). — J dxC{x)U{x) -f J dxdyD{x,y)U{x)U{y). Ir ' must be zero for any symmetric potential. 
Substituting U{x) = US{x — xo) we find that C(xo), D{xo,xo) = 0. Substituting U{x) = US{x — xi) + US{x — X2) we 
see that D{xi,X2) + D{x2,xi) = 0. Hence, ' = for any U{x). 

We first consider the 'non-trivial' ratchet effect and then check what changes after the voltage dependence of the 
injected charge density is taken into account. Let us begin with a qualitative explanation before we make a rigorous 
calculation. The origin of the ratchet current can be understood from a simplified Hartree-Fock picture. In this 
approximation, electrons are backscattered off a combined potential U{x) = U{x) + W{x), where W{x) is a self- 
consistent electrostatic potential created by the average local charge density. To obtain W(x) we use the following 
approximation in the last term of Eq. JTJ: p{x)p{y) « (7/>^(a;)V'fi(a;) -I- 7/;^(a;)V'L(a;))(V'fl(y)V'-R(y) + V'l + 
[{p{x))ifi^{y)'ipL{y) + {p{y))^R{x)ip + h.c] + const. Thus, the relation between W and p is linear. The combined 
potential U{x) is different for the opposite voltage signs. 

In the model the electron interaction is short-ranged due to the screening gates, and hence, the relation between 
the potential W{x) and the electron density p{x) is local, W{x) ~ p{x). The simplest choice of U{x) is a two-impurity 
asymmetric potential U{x) — UiS{x + a/2) + U2S{x — a/2). The charge density profile in the presence of a two- 
impurity potential and the voltage drop V is sketched in Fig. ^ Depending on the voltage sign, the charge density 
decreases or grows as a function of the coordinate x. So does the electrostatic potential W{x). Hence, U{x) is different 
for the opposite voltage signs. The density is essentially independent of the coordinate between the impurities 0|, as 
well as on the left and on the right of the impurities, since no backscattering occurs in those regions. The charge density 
and the electrostatic potential drop at the positions of the impurities. The magnitude of a drop is proportional to the 
electric charge backscattered off the impurity. Indeed, if the incident charge densities of the electrons approaching 
the impurity from the left and from the right are pj/ and p*^ , and the backscattered charge densities are pj^ and p~^ 
then the density drop across the impurity Ap = {pj/ + p*^ + p*^ — p~^) — (pj^ -I- p~j^ + p~j^ — p*£) = 2(p^ — p~j^) ^ /bs, 
where /bs is the current backscattered off the impurity. Thus W{x) ^ U — U ~ /bs- From Ref. 0| we know that for 
a weak potential U 

hs-\U2kA'\V\^'-'signV/El^, (2) 

where U2kF ^ kp J dx e'xjp{2ikpx)U (x) , kp is proportional to the mean electron density, and the dimensionless constant 
g characterizes the interaction strength, g — I for non-interacting electrons (in which case W{x) = 0). 

Now we can substitute the renormalized potential U = U + W for U in Eq. The Fourier component W2kF 

is different for the opposite voltage signs. Hence, we obtain the asymmetric part of the I — V characteristics ~ 
eU^\eV\'^^~'^ / (hE^). The ratchet effect is strongest for 17 ^ when the ratchet current grows as the voltage decreases. 

The above Hartree-Fock argument provides a qualitatively correct picture at small g but underestimates fluctuations 
in Luttinger liquids. As shown below, the ratchet current growth at small voltages differs from our estimate: Ir ^ 
f/3|y|6g-2^ Ep :$> V > V* {U E~^)^/ ^^s) ^ 5 < 1. We will see that the growth terminates at F = V* . At 
such voltage Ir{y*)/I{V*) - [{y*)^a+i / e]^]/[v*] ~ {V*/Ep)^9 ~ 1 as g < 1. Fluctuations are less important in 
many-channel systems and the Hartree-Fock picture gives exact results for some two-channel systems and for Fermi 
liquids 0. 

We use the bosonization technique to calculate the ratchet current. After an appropriate rescaling of the time 
variable, the system can be described by the action 

S= /"didx|^Pt$)2-(c),$)2]_5(a:)^2C72„fc,cos(n^$ + a„)|, (3) 

^ n>l ^ 

where the bosonic field $ is related to the charge density as p = e{y/gdx'^ + 2kp)/{2TT), and U2nkF exp(ia„) are of the 
order of the Fourier components of the asymmetric potential, kp J exp{2inkpx)U {x)dx. We assume that the charge 
density ~ kp is independent of the voltage. The operator cos(ny^$ -I- a„) describes scattering events involving n 
electrons. We assume that ai = 0. Indeed, we can always set ai = by a constant shift of the bosonic field $. For 
a general asymmetric potential, a„ with n > 1 remain non-zero after this shift. On the other hand, for a symmetric 
potential U{x) = U{—x) all a„ — 0. In most problems it is sufficient to keep only the n — 1 term. The n — 2 
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contribution is relevant in the theory of resonant transmission in Luttinger hquids |13| . This term is also important 
for the ratchet effect. 

We use the standard model [l8lll9ll2oj for Fermi-liquid leads adiabatically connected to the wire. We assume that 
the action Q is applicable for \x\ < L only. At large |a;| the interaction strength K{x — y), Eq. is zero. This 
model can be interpreted as a quantum wire with electron interaction completely screened by the gates near its ends. 
Electric fields of external charges are assumed to be screened in all parts of the wire. A simple modification of this 
model describes electrically neutral leads . All results coincide for our set-up and the model [2^ . 

The current injected from the non-interacting ID regions is given by the Landauer formula Iq = e^V/h [l^ . Indeed, 
lcft-/right-movers entering the non-interacting region from the central part of the wire cannot affect the current of 
right-/left-movers in the non-interacting region. Hence, the current of right-/left-moving particles in the left/right 
non-interacting region is determined by the chemical potential of the left/right reservoir. The total current is the sum 
0,|23| of the injected current and the current backscattered off the asymmetric potential: I = Iq + I^s- Only /bs 
contributes to the ratchet effect. To find the backscattered current we employ the Keldysh formalism ^2)\. We assume 
that at t — —00 there is no backscattering in the Hamiltonian, and then the backscattering is gradually turned on. 
Thus, at the initial moment of time the numbers Nj^ and Nj^ of left- and right-moving electrons conserve separately. 
Hence, at t = —00 the system can be described by a partition function with two chemical potentials (Xr = Ep and 
/iL — Ep + eV conjugated with the particle numbers Nj^ and N^. This initial state determines bare Keldysh Green 
functions. 

We will consider only zero temperature. It is convenient to switch jl9l | to the interaction representation H — > 
H — fj,RNii~ (ilNl- This transformation induces time dependence in the electron creation and annihilation operators. 
As a result, J^n '2U2ri,k^ cos(n^/^^ + a„) in the action should be modified as yi„ '^ U-inkip cos(ny^$ + a-n + 'nA(t)), 
where A{t) = eVt/h ^^E^ 

The backscattered current operator equals |l3l |20j 

hsit) = dNL/dt = i[H, NL]/h - -6S/SAit), (4) 

where we omit dimensional factors such as e, h, and vp for brevity. We need to calculate 

(/bs(i = 0)) = (0|5(-oo; 0)/bs(0)5(0; -oo)|0), (5) 

where |0) denotes the initial state and S is the evolution operator. In the weak impurity case this can be done with the 
perturbation theory in (l2nkF using the bare Green function {0\^{ti,xi = 0)^{t2,X2 = 0)|0) = —2\n{S + i[ti — t2]), 
where d is an infinitesimal positive constant. 

If all ak — then the ratchet current is zero. Indeed, at ~ the action (O is invariant under the transformation 
$ — $, V —V while the current operator |0J changes its sign. As discussed above, for an asymmetric potential we 
expect a2 7^ 0. Then a ratchet current emerges in the order U^,,^Uik„. Before the calculation of 7^. let us determine 
its voltage dependence with a heuristic argument similar to Ref. [l3|. As one changes the energy scale E, the 
backscattering amplitudes U2nkp in the action Q scale as ij2nkF{E) ~ U2nkFE'^^s~^ ^M- This renormalization stops 
at the energy scale V. Assuming that a scattering matrix approach could be applied for an estimation of the current, 
we write 4.(1^) - VR^s{V), where Res{E) = E constt/|„^^^ (7;) +J2constij2nkAE)U2mkAE)U2ikAE) + ... is an 
effective reflection coefficient. Quadratic terms do not contribute to the ratchet current. The leading contribution 
emerges in the order tj^f^^UAkp- One gets Ir ^ ^f^lfc^ (^)^4'^f (^) ~ V^^~'^ . Below we obtain the same result 
rigorously from Eqs. H4I5() . 

Expanding Eq. (0) to the order U^^kF gives 

Ir = 2sina2t/2fc^^4fcf ( / dti f dt2Cos{Vti - 2Vt2)P{tiMM ~ h) 

dti dt2Cos{Vti~2Vt2)P{ti,t2,h-t2)- dh dt2Cos{Vti-2Vt2)P{-h,t2-h,t2) 

-2/ dti I dt2Cos{Vti+Vt2)P{t2-ti,tiM)+ I dh dt2Cos{Vti + Vt2)P{t2-ti,-ti,t2)\ +c{6: 

where P{t, s, q) = {6 — itY^{5 — is)^'^^{6 — iq)^'^^ . Dimensional analysis shows that Ir ^ V^^^^ in agreement with our 
previous estimate. It is convenient to change variables in the integrals with cos(yii 2Vt2) as ti = t2 — ti, T2 = t2. 
Then after tedious but straightforward manipulations Eq. ^ can be represented as 

f pOG poo ^ 

Ir = 2sma2Uiki.U4kF \ dtdr cos[V{t - t)]P{-{t + t), -t, -t) - dtdr cos[V{t + t)]P{T - t,T,t)[ . (7) 
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The first integral in ^ is zero as seen from the location of the branching points of the function P. The second integral 
yields 

= - .inaA^U,,, cos(7rg)^^^^^?|^^^|±^ (8) 

This expression becomes at g = 1/2. We also get a zero ratchet current for non-interacting electrons, g = 1, 
because the Hamiltonian is quadratic in Fermi-operators in the non-interacting case and hence no operators which 
backscatter more than one electron can appear, U4kp = 0. 

At small g the ratchet current Q is proportional to a negative power of the voltage. This means an unusual 
behavior: the dc response to an ac voltage grows as the ac voltage decreases. 

So far we ignored the voltage dependence of the injected charge density. At g ^ 1, Eq. © gives the main 
contribution to the ratchet current only for eV < \/UEp. For g close to 1 the result ^ is always exceeded by another 
contribution. This contribution emerges in the second order in U and is related to the voltage dependence of the 
injected charge density. The density is proportional to kp which enters the expression for U2kF iii Eq. At small 
V <C Ep the correction to U2kF is a linear function of V . The substitution of this correction into Eq. Q gives 
an additional ratchet current 

, , . eC/2, (eV fs 

For g > 1/3 and V > V* the contribution ^ always exceeds At 5 < 1/3 the current ^ is greater than 
© above a threshold voltage that depends on U and g. As we already discussed, Ir ^ is comparable with the 
total current I{V) ^ e^V/h at small g near the border of the perturbatively accessible region UV^~^ / Ep < 1. On 
the other hand, Eq. provides only a small correction to the total current for any g. Still a repulsive interaction 
of any strength enhances the ratchet effect as seen from the comparison of the current ^ for g < 1 and for the 
non-interacting case 5=1. 

What happens beyond the perturbative region when V <V* ^ As the energy scale decreases the effective 

impurity strength grows. Hence, we need to consider a strong U > Ep limit. In this limit we have a weak tunneling 
between the left and right halves of the wire. The current I{V) ~ t'^V^'^^^'''^ / Ep^^ , where t is the tunneling amplitude 
[l3| . Inserting the voltage dependence of the tunneling amplitude in the expression above we estimate Ir (V) ^ l/^/s _ 

A single impurity model Q can be used only when the potential U(x) is confined in a small space region of 
size a < av ^ Tivp/{eV). If the potential changes slowly at the scales x > ay ^ 1/kp it cannot backscatter 
electrons since backscattering involves high momentum transfers, AA: > kp. Interesting interference effects are 
possible for a two-impurity potential Uid{x) + U26{x — a) and other U{x) which significantly change at the scale 
l/kp but are non-zero in a region of size a ay- In the two-impurity case the current oscillates as a function of the 
voltage bias For t/i, U2 Ep, I - e^V/h ~ [Ul + C/| -I- 2U1U2 cos{2kpa)H{geVa/[nvp])]\V\^3-isyg^v , where 
H{x) = 0rr(25)Jg_i/2(a;)/[r(g)(2x)f~i/2] Jg^i/2{x) is the Bessel function of the first kind [13. The main 
contribution to the ratchet current at a ~ ay comes from the shift oi kp due to the change of the electrochemical 
potential of the left reservoir by eV . From the minimum of the quadratic part of the bosonized Hamiltonian one finds 
the charge density shift This gives kp — k^p' + g"^ eV / {2hv p) . After the substitution to the expression for the 

total current / we find 

Ir{V) - C/iC/2sin(2fc^°^a)|T/pf-isin(g2e|l/|a/[?iwF])i?(5eya/[?iui.]) (10) 

Thus, Ir{y) oscillates. Notice that for V ^ V* ^ Ep, a ^ ay the ratchet current (|10|l is of the order of the total 
current ^ e^V/h. 

In conclusion, we have found the ratchet current for strong and weak asymmetric potentials. It exhibits a set of 
universal power dependencies on the voltage and can grow as the voltage decreases. This work was supported by the 
US DOE Office of Science under contract No. W31-109-ENG-38. 
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FIG. 1: Density profiles averaged over the period of Priedel oscillations for a potential with Ui < U2- The averaged densities 
show drops at the impurity positions. 



